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We investigate non-perturbative results of inviscid forced Burgers equation supplemented to continu- 
ity equation in three-dimensions. The exact two-point correlation function of density is calculated 
in three-dimensions. The two-point correlator < /9(xi)p(x2) > behaves as |xi — X2|~°^ and in 
the universal region as = 7/2 while in the non-universal region aa = 3. In the non-universal re- 
gion we drive a Kramers-Moyal equation governing the evolution of the probability density function 
(PDF) of longitudinal velocity increments for three dimensional Burgers turbulence. In this region 
we prove Yakhot's conjecture [Phys. Rev. E 57, 1737 (1998)] for the equation of PDF for three 
dimensional Burgers turbulence. We also derive the intermittency exponents for the longitudinal 
structure functions and show that in the inertial regime one point Urma enters in the PDF of velocity 
difference. 

PACS numbers, 47.27.Gs and 47.40.Ki 



I. - INTRODUCTION 

Recently, tremendous activity has started on the non- 
perturbative understanding of turbulence [1-22]. A sta- 
tistical theory of turbulence has been put forward by Kol- 
mogorov |2^, and further developed by others [p^-p6[. 
The approach is to model turbulence using stochastic 
partial differential equations. The simplest approach 
to turbulence is the Kolmogorov's dimensional analysis, 
which leads to the celebrated k~^/^ law for the energy 
spectrum. This is obtained by decreeing that the en- 
ergy spectrum depends neither on the wavenumber where 
most of the energy resides, nor on the wavenumber of vis- 
cous dissipation. Kolmogorov conjectured that the scal- 
ing exponents are universal, independent of the statistics 
of large-scale fluctuation and the mechanism of the vis- 
cous damping, when the Reynolds number is sufficiently 
large. In fact the idea of universality is based on the no- 
tion of the " inertial subrange" . By inertial subrange we 
mean that for very large values of the Reynolds number 
there is a wide separation between the scale energy input 
L and the typical viscous dissipation scale 77 at which vis- 
cous friction become important and the energy is turned 
into heat. 

In turbulent fluid flows, the physical quantities of inter- 
est, such as the velocity and density fields, display highly 
irregular fluctuations both in space and time. To study 
such processes, using a statistical approach is most nat- 
ural. The statistics of any fluctuating quantity are de- 
scribed by its probability density function (PDF). One 
would, therefore, hope to be able to calculate the PDF's 
of the turbulent quantities directly from the equations 
of motion, i.e. for example to calculate the PDF of ve- 
locity fluctuations from Navicr-Stokes equation. How- 



ever this task is highly nontrivial and as of today has 
not been accomplished not even for the relatively simpler 
problems. The striking feature of developed turbulence 
is its intermittent spatial and temporal behavior. The 
structures that arise in a random flow manifest them- 
selves as high peaks at random places and at random 
times. The intervals between them are characterized by 
a low intensity and a large size. Rare high peaks are 
responsible for PDF tails while the regions of low inten- 
sity contribute PDF near zero. Analytically the follow- 
ing properties manifest themselves in the nonlinear ba- 
havior of the scaling exponents in the velocity difference 
structure functions (transvers and longtudinal) 5„(r) i.e. 
< {u{x + r) — u(a;))" >. It is known that S'„(r) behaves 
as ~ r^" , where C„ is a non-trivial function of n. This 
behaviour leads to stongly non-gaussian (intermittent) 
probability distribution functions of velocity difference 

On the other hand Burgers equation, which describes 
the potential flow of the fluid without pressure, provides 
a wonderful laboratory for testing new ideas and tech- 
niques in view of the study of fully developed turbulence 
in the Navier-Stokes turbulence. These are two cases of 
non-linear stochastic equations which share in the same 
structure of the non-linearity. The important differences 
come from the nature of interaction (i.e. locality and non- 
locality) and the large scale structure. In the case of in- 
compressible turbulence the interaction is non-local that 
is the very existence of the pressure causes the far regions 
of the ffow to be coupled with together ]jll| and leads to 
effective energy redistribution between components of the 
velocity field, while in the case of Burgers equation the 
pressure effects are absent. The non-linear term (v • V)v 
in the Navier-Stokes and Burgers equations tends to form 
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shocks and enhance intermittency while the additional 
pressure gradient term in N-S equation depresses inter- 
mittency. This senario has been given addressed by V. 
Yakhot [Tl| ] . for constructing an approxiamte dynamical 
theory for N-S equation. In multi-dimensional Burgers 
turbulence the presence of large scale structures (shocks) 
forming a d-dimensional forth-like pattern is responsible 
for extreme case of intermittency causing the saturation 
of the intermittency exponent to Cn — 1- Large scale 
structures in "true" turbulence are similarly thought to 
be the origin of the experimentally observed intermit- 
tency, which is however much milder. The deep reason 
of this difference is thought to be probably related to the 
"dimension" of the large scale singularities, which is d — 1 
in the d-dimensional Burgers turbulence and only 1 for 
votex lines in hydrodynamical turbulence. 

The problem of forced and un-forcerd Burgers turbu- 
lence 28-30] has been attacked recently by various meth- 
ods P ,3~|6|,P|~^ ■ The Burgers equation descibes a vari- 
ety of non-linear wave phenomena arising in the theory 
of wave propagation, acoustics, plasma physics, surface 
growth , charge density waves, dynamics of the vortex 
lines in the high superconductor, dislocations in the 
disordered solids and formation of large-scale structures 
in the univers e |]30| - |3^ . According to the recent theoret- 
ical |4p|, pT| , [r7| , [l8| and numerical works [Q, it is known 
that the PDF for the velocity difference behaves differ- 
ently in universal and non-universal regions. In the uni- 
versal region i.e. the interval | Am| << Urms and r « L, 
the PDF can be represented in the universal scaling form 
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(1) 



where F{x) is a normalizable fuction and the exponent z 
is related to the exponent of random force correlation 77 
as z = ^4-^ ■ In the region for x — -^^^ > > 1 the universal 
scaling function F{x) is given by the expression F{x) ^ 
exp(— aa;'^), where a is some constant in one-dimension 
and it depends on the cosine of angle between the vectors 
Au and r in the higher dimensions. On the other hands 
the PDF in the interval |Au| >> Urms behaves as: 



P(Au,r)=rG(-^) 



(2) 



which depends on the single-point Urms and therefore is 
not a universal function. 

In this paper we consider the inviscid forced Burg- 
ers equation in 3-dimensions supplmented to continu- 
ity equation. We find the exact two-point correlation 
functions of density in three-dimensions and show that < 
p(xi)p(x2) > behaves as |xi — X2|~"^, where as = 7/2 
in the universal region i.e. \u{x) — u{x')\ « Urms where 
rms value of velocity fluctuations. In the 
non-universal region i.e. \u{x) — u{x')\ » Urms we drive 
the exact Kramers-Moyal equation governing on the evo- 



lution of the probability density function (PDF) of longi- 
tudinal velocity increments for three dimensional Burgers 
turbulence. In this region we prove the Yakhot conjec- 
ture for equation of the PDF in multi-dimensional 
Burgers turbulence and prove the existence of the A-term 
and B-term in the equation for the PDF of longitudinal 
velocity difference. We also derive the intermittency ex- 
ponents for the longitudinal structure functions and show 
that how the Urms enters in the PDF of the velocity dif- 
ference and breaks the Galilean invariance. Furthermore 
we refer to the zero viscosity limit of 3D forced Burgers 
turbulence. We claim that at least when the simplest 
terms satisfying the basic symmetries of the problem are 
introduced for closing the disspiation term the require- 
ments of the homogeniety and isotropy would force the 
cocfficcient of the proposed terms to vanish. So within 
this approximation the viscous term in the vanishing vis- 
cosity limit do not contribute in the equation for the 
PDF's of velocity difference in the inertial range. 

The paper is organized as follows: In section 2 we de- 
rive the density-density correlators exponent in invicsid 
Burgers turbulence and the generating function for lon- 
gitudinal velocity difference in the universal part of the 
PDF . In section 3 we determine the small scale statisics 
of longitudinal velocity difference and derive the exact 
equation for the PDF's of longitudinal velocity difference 
in the non-universal part of the PDF and derive the ex- 
act value of intermittency exponents in longtidinal struc- 
ture functions. We show that the intermittency problem 
in these systems can be investigated non-perturbatively 
and derive the deformation of PDF's in length scales. It 
is shown that deformation of PDF from large to small 
scales are completely described by Kramers-Moyal equa- 
tion. 



II. - PDF'S OF 3D BURGERS TURBULENCE IN 
THE UNIVERSAL REGION 

Our starting point is the 3D-Burgers equation supple- 
mented to continuity equation: 



Ut (u • V)u = vV'^u + f (x, t) 



Pt + da[pUa) = 



(3) 



(4) 



for the Eulerian velocity u(x, t) and viscosity v and den- 
sity p, in 3-dimensions. The force f(x, i) is the exter- 
nal stirring force, which injects energy into the system 
on a length scale L. More specifically one can take, for 
instance a Gaussian distributed random force, which is 
identified by its two moments: 



< /^(x,<) > =0 
< /^(x,i)/,(x',t') > =fc(0)(5(i-i')V(x- 



(5) 
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where fx^v = x,y,z. The correlation function k^y{r) 
is normaUzed to unity at the origin and decays rapidly 
enough where r becomes larger or equal to integral scale 
L. The quantity fc(0) measures the energy injected 
into the turbulent fluid per unit time and unit volume, 
f (x, t) provide also the energy flux in the k — th shell as 
Ilfc = n(r ~ J^^j^ < |f(k)p >, where r belongs 

to the inertial range. 

Eqs.(3) and (4) exhibit special type of nonlinear inter- 
actions. It is hidden in the nonlinear term (u • V)u. The 
advective term couples any given scale of motion to all 
the large scales and the large scale contain most of the en- 
ergy of the flows. This means that the large-scale fluctu- 
ation of turbulence production in the energy-containing 
range couple to the small-scale dynamics of turbulence 
flow. In other words, the details of the large-scale tur- 
bulence production mechanism are important, leading to 
the non-univesality of probability distribution function 
(PDF) of velocity difference. However in the case that 
\u{x) — u{x')\ << Urms it is belicvcd that the PDF for the 
velocity difference is not depended to Urms and therefore 
the one-point Urms does not appear in the velocity differ- 
ence PDF. This region is known as the Galilean invariant 
(GI) region. 

Now let us consider the invisid 3D forced Burgers tur- 
bulence that is we intend to neglect from the viscosity 
contribution in eq.(3). The problem is to understand the 
statistical properties of velocity and density fields which 
are the solutions of eq.(3) and (4). We consider the fol- 
lowing two-point generating function: 

^2(Al, A2,Xi,X2) = 

< p(xi)p(x2) exp(Ai.u(xi) + A2.u(x2)) > (6) 

where the symbol < • • • > means an average over various 
realizations of the random force. 

To derive an equation for F2, we write the eq.(3) 
and (4) in two points xi and X2 for ui,U2,ut, and p(x) 
and multiply the equations in p(x2), Ai2;p(xi)p(x2), 
Ai2p(xi)p(x2) and p(xi), A2xP(xi)p(x2), • • •, and 
A2zp(xi)p(x2), respectively. We add the equations and 
multiply the result by exp(Ai.u(xi) -I- A2.u(x2)) and 
make average with respect to external random force, so 
we find: 



dtF. 



E 



d 



{?'— 1.2}/j-— a:,y,j3 

E Ai,^A-,,i,A:^t.(xi - Xj)F2 = (7) 

— l,2}/i,i^— a;,y,2 

where we have used the Novikov's theorem. The above 
equation is first driven by Polyakov It is noted that 
the advection contributions are accurately accounted for 
in this equation. Also we suppose that kfn, has the fol- 
lowing form: 



V(xi - Xj) = fc(0)[l - ^^i^^^^,, 
_ (xj ~ x:j)p(xi — Xj)^ 

with k{0),L — 1. We change the variables as: x± 



(8) 



Ai-A 



- and and con- 



xi ± X2, A+ = Ai + A2 and A_ = 
sider the spherical coordinates, so that X- : (r, 9, (f) and 
A- : {^,9',(p'). Hence we find that the F2 satisfies the 
following closed equation for homogeneous and isotropic 
Burgers turbulence: 



[sdrdf. 
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-ds 



drds 



(9) 



where s — cos 7 = cos 6* cos 6*' + sin 6' sin 6*' cos((y5 — </?'). 

In homogeneous and isotropic Burgers turbulence, 
with stirring correlation as A:(r) ^ 1 — r'' (where in the 
our case i.e. eq.(8) we have ry = 2), we consider the uni- 
versal scale-invariant solution of eq.(9) in the following 
form: 



^2(m, r) ^ g{r)F{nr ) g{r) ^ r~ 



(10) 



Substituting the following form for the generating func- 
tion fixes the exponent 5 as S = ^^ii (In our case using 
eq.(8) we find 6 = 1). Invoking to the scaling invariance 
of the inviscid Burgers equation and continuity equation, 
we assume the existence of the density-density correla- 
tors with the scaling form of introduced in (10). as is 
the exponent of two point correlation and the two point 
correlation of the density can be found by the generating 
function in the limit of /i — > 0. Therefore it is neces- 
sary to find such a solution for F{^r^) which tends to a 
constant in the limit of 0. Proceeding further we 
focus our attention to the longitudinal velocity compo- 
nents, i.e. s = 1, and hence accept the scaling ansatz 
F{iir,s) = F{fj,rs). The proposed form of the argu- 
ments will garantee that S'„(r, s) ^ s"S'„(r) when n < 1. 
Rewriting eq.(9) in terms of the variable z = firs, the 
following equation is obtained in the limit of s ^ 1: 



zd'^F{z) + (3 - a3)d,F{z) ^ 3z^ = 



(11) 



It is interesting that the above equation was first de- 
rived by Polyakov 1^ for the problem of one dimensional 
Burgers equation in the inviscid limit. In that work the 
effect of the viscose term is found in the limit of 
and r ^ L by appeal to the self-consistent conjecture of 
operator product expansion. It is found that there are 
only two terms generated by the viscose term which are 
consistent with the symmetries of the problem. These two 
anomaly terms will modify the master equation governing 
over the generating function in a way such that a positive 
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, finite and renormalizable PDF can be found [|. A sim- 
ple analogy between eq.(lO) and Polyakov's result will re- 
veal that the coefiicient of b- anomaly is replaced with the 
scaling exponent of density-density correlation through a 
simple linear relation. In the problem of one-dimensional 
Burgers equation in the zero viscosity limit the presence 
of the b- anomaly generated by viscosity term is crucial 
for finding a positive PDF for velocity increments in the 
universal regime and the requirement of positivity will 
fix the value of anomaly coefficient ||^. Boldyrev p^ ] 
shows that one could find a family of solutions for dif- 
ferent values of the b-anomaly coefficient if one relaxed 
the homogeniety condition for the universal part of the 
PDF. The value of this coefficient is related to the alge- 
braic decay of the left tail of PDF in the universal regime. 
Determination of that decay exponent has been a contro- 
versial subject for which other methods have been devel- 
oped and among them recent rigorous methods should 
be mentioned within which it is fixed to 7/2 |l^Jl^. The 
interesting point is that our calculations in three dimen- 
sions show that when density fluctuations are taken into 
account there is no more any need to the viscosity term 
for obtaining a positive PDF and even in the inviscid 
problem it is possible to find a positive solution for the 
longitudinal velocity increment PDF. It is easy to show 
that the requirement of the positivity on the PDF will fix 
the density-density scaling exponent to — 7/2. There 
are some strong constraints which would be considered 
for any non-perturbative dynamical theory constructed 
for describing the intermittent behaviour of the structure 
functions in the homogenious and isotropic turbulence 
fiows. The most vital constraints are the homogeniety 
and isotropy which dictate the average of the velocity in- 
crement to be vanished. Now it is straight forward to 
show < u >s= 0, where u = u(xi) — u(x2) [|9|. 

We have also extended our analysis to two-dimensional 
Burgers turbulence. We define the variables as: x± = 
xi ± X2, A+ — Xi + A2 and A_ — ^^^^^ ^"^^ consider 
the polar coordinates, so that x- : (r, 6) and A_ : (/i, </?). 
Now one can show that for isotropic and homogeneous 
Burgers turbulence, two-point generating function F2 
satisfies the following equation: 

r/i r/i /i 

+ ^—^d^ds-r^fi^il + 2s^)]F2^0 (12) 
r 

where s = cos{6 — ip). Similar to three-dimensions, 
we propose the scale invariant solution for F2 as F2 = 
g{r)F{z, s) where g{r) = r~"^ and z = rfi. In the limit 
of s = 1 once again it would be natural to seek the scaling 
solutions as F{r^, s) = F{rsfi). It will be straightforward 
to find that the eq.(ll) will again governs over F2{firs) 
and the only difference originates from the coefficient of 
first derivative with respect to z = ^rs which in this case 



is 2 — a2- Using the positivity and normalizability con- 
dition of PDF we find F{iJ,rs) = exp(2:''') and show that 
a2 — 5/2 and 7 = 3/2. The exact values of the density- 
density exponents i.e. 02 and 03 are the main results of 
this part of our work which is derived self consistently. 

Since the velocity difference PDF is the Laplace trans- 
formed of F, one can readily deduce the right tail of the 
PDF as ^ exp(-(|^)3) (for s = 1) in two and three di- 
mensions and in the limit ^ —i- +00. This tail has been 
confirmed by several other approaches [pp|Jl5[|. Left tail 
of the PDF is sensitive to the scaling exponent of the 
density-density correlator and is given by (^„)(q3-i) when 

— 00. At this stage we can not derive the PDF of 
density fluctuation. The reason is that wc have only two- 
point correlation functions of the density flcld while all 
of moments of density field are needed for this purpose. 
In the one dimensional decaying Burgers turbulence it 
has been claimed that the density PDF would have some 
power law tail |l2| but this problem is open for the forced 
case. 

III. - PDF'S OF 3D BURGERS TURBULENCE IN 
NON-UNIVERSAL REGION 

In this section we consider the 3-dimensional Burg- 
ers turbulence in the non-universal region i.e. \u{x) — 
u{x')\ >> Urms, and derive the PDF of the longitdi- 
nal velocity difference and therefore the exponent of ve- 
locity structure functions. As mentioned in the intro- 
duction recent works [^|J^, [ll| , [l7| , |l8| indicate that in the 
non-universal region the PDF of velocity difference de- 
pends on the one-point Urms and therefore is not univer- 
sal which is meant to be sensitive on the details of large 
scale forcing. This problem is known as the break down 
of Galilean invariance in the non-universal region. The 
force free Burgers equation is invariant under space-time 
translation, parity and scaling transformation. Also it is 
invariant under Galilean transformation, x x + Vt and 
V — > v + V, where V is the constant velocity of the moving 
frame. Both boundary conditions and forcing can violate 
some or all of symmetries of force free Burgers equation. 
However it is, usually assumed that in the high Reynolds 
number flow all symmetries of the Burgers equation are 
restored in the limit r — > and r >> 77, where rj is the 
dissipation scale where the viscous effects become im- 
portant. This means that in this limit the root-mean 
square velocity fluctuations Urms = V < v-^ > which is 
not invariant under the constant shift V, cannot enter 
the relations describing moments of velocity difference. 
Therefore the effective equations for the inertial-range 
velocity correlation functions must have the symmetries 
of the original Burgers equations. For many years this 
assumption was the basis of turbulence theories. But 
based on the recent understanding of turbulence, some 
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of the constraints on the allowed turbulence theories can 
be relaxed §,|ll|,0[l8|. 

According to recent Yakhot modeling of Burgers and 
N-S turbulence this symmetry breaks in a hard way i.e. 
the Urms is entered explicitly in the equation of PDF 
for velocity difference. In the following we aim to show 
how this symmetry breaks in the sense that the one-point 
Urrns cnters in the argument of the PDF in non-universal 
region. Also since we are interested in the scaling of 
the longitudinal structure function Sq =< {u{x + r) — 
u{x))'^ >=< u'' >, where u{x) is the x-component of the 
three-dimensional velocity field and r is the displacement 
in the direction of the x-axis and the probability density 
P(u, r) for homogeneous and isotropic turbulence. In 
the non-uiversal region using the eq.(ll) one can observe 
that the PDF for the velocity difference in the invicsid 
Burgers turbulence in two and three dimensions satisfy 
the following closed equation: 



—^duuP — sduudrP ^- 

r r 



d-2 



-dsP 



ad 



{1 ~ s^)dsP + (l - s^)dsdrP 



duudsP 



r2(l 



2s^)dlP = 



(13) 



where s equal to cos{6 — ip) in two-dimensions and 
s — cos 7 = cos 6* cos 6*' -I- sin0sin6''cos((p — cp') in three- 
dimensions. For determining the statistics of the small- 
scale and the behaviour of PDF in the inertial range, i.e. 
r] << r << L, we can ignore the forcing term in eq.(19). 
The reason for neglecting the force term is that the forc- 
ing contribution is in the order of r'^ /L"^ so in the inertial 
range we can safely drop the corresponding term. How- 
ever we need to consider the effect of forcing by matching 
the PDF in the inertial range with the PDF in the in- 
tegral scale. Since we are interested in equation of the 
longitudinal PDF we have to consider the above equa- 
tion in the limit of s — > 1 . This limiting is not trivial and 
needs to be considered more carefully. The contributions 
of different terms of the PDF equation in the limiting 
when s — !■ 1 is determined by the corresponding equation 
of the structure function. 

Assuming that the all of the moments of velocity dif- 
ference exist, the structure functions Sn for given angle 
7 ( or s = cos(7)) satisfies the following closed equation: 

[sn -f (1 — s'^)ds\rdrSn ~ nsadSn 
-s{l - 5„ + [d-2 + s^)dsSn 
-ad{l - s^)dsS.n + n{l - S^)dsSn 



+r^n{n - l)(n - 2)(1 + 2s^)Sr,-z = 



(14) 



The forcing contribution to the above equation is the last 
term i.e. r'^n(n— l)(n — 2)(l + 2s^)5„_3 and this term has 
not any contribution in the exponent of structure func- 
tion. However the amplitude of the structure functions 



do depend on the details of forcing. We mean that the 
scaling ( multifractal ) exponent of the structure func- 
tion is not related to the forcing term and is related to 
the structure of non-linearity and the transverse contri- 
butions in the Burgers equation. For solving the eq.(20) 
one can seperate the angular and scale dependent parts 
of S'„(r, s). The calculations give rise to the result that 
for s \ the structure functions have the following form 



5„(r, s)^s"5„(r), 



(15) 



where S'„(r) =< (u(a; + r) — u(a;))" >=< u" >. Factoriz- 
ing the angle and scale dependences is known also for the 
N-S turbulence too . The proposed form for struc- 

ture functions dictates that in the limit when s — » 1 the 
probability distribution of velocity increments satisfies 
the scaling form, P(r, u, s) — > ■^P{r,u/s) in a sufficient 
way. It is clear that when s — > 1, P{r,u) = P{—r,—u) 
and it satisfies the following equation in d=2 and 3- 
dimensions. 

ou or r ou ou'^ 

where P is the longtudinal velocity difference PDF, and 
B approaches to zero as 0(1 — s^) and giving address to 
(■3 = 1 would give A=l. The above equation is correct 
for three and two dimensions and it does not depend on 
the dimensionality. This form of equation for the PDF 
has been conjectured recently by Yakkot for the multi- 
dimensional Burgers turbulence [|ll[ . It is noted that the 
forcing contribution in the above equation is 'ir^d^P and 
it is irrelevant in the small scale r — + 0. We will deal with 
the forcing contribution as a boundary condition of PDF 
in the large scales and this term is responsible for the 
breakdown of the Galilean invariance . Also in the non- 
universal region we find the density-density exponent to 
he ad = d in the two and three-dimensions. Now it is easy 
to see that the eq.(8) can be written as a scale-ordered 
exponential: 

I dr' LKM{n,r') 

P{u,r)^T{e+° P{u,ro)) 

where Lkm can be obtained by computing the inverse 
operator |34|. Using the properties of scale-ordered ex- 
ponentials the conditional probability density will satisfy 
the Chapman-Kolmogorov equation. Equivalently we de- 
rive that the probability density and as a result the con- 
ditional probability density of longitudinal velocity incre- 
ments satisfies a Kramers-Moyal evolution equation: 



dP 
dr 



n=l 



hir^{D(-\r,u)P) (17) 



Where D^"'>{r,u) = ^u" @. We find that the co- 
efficients (T„ depend on A and B through the relation 
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an = (~l)" (B+l)(B+2)(1?+2)...(B+n) - ^^^^ ^q^^^" 

tion (i.e. eq.(23)) obviously governs over the condi- 
tional PDF too but with another boundary condition, 
i.e. P{u,r\u' ,r) — S{u — u'). For a simple case we pro- 
ceed to find velocity difference PDF P{u, r) just by using 
the same line of reasoning and the well known Bayesian 
rule. So 



P(M,r)= P{u,r\u',L)P{u',L)du' 



(18) 



where P{u',L) is assumed a Gaussian in the integral 
scale. Since 

/ dr'LKM{U,r') 

P{U,r\u',r) =T{e_;" S{u - u')) (19) 

The proposed KM operator has an important prop- 
erty,i.e. 



Cm"" 



(20) 



where Cm is the scaling exponent of the longitudinal ve- 
locity difference 5m. Plugging the scale ordered form of 
the conditional PDF in the eq.(24), we will get: 



P{u,r) = iy)^KM(^)p(u,L) 



(21) 



Expanding the assumed gaussian form of P(w, L) in 
terms of u and using (26) ends with, 



exp[ln(f )C2m] , U 2m 



m=0 



^ rras 



7 2^ 



" )2m 



^ ^ rras 

where Cj — 1- We should stress that the same structure 
is tractable for the other PDF's in the integral scale other 
than the simple Gaussian form. This result could be de- 
rived from eq.(22) by direct calculations even without 
consulting the KM form of the evolution operator too. 
This form verifies the proposed form of the PDF in the 
inner scales where r\ « r « L and resonates with nu- 
merical observation in the one dimensional Burgers tur- 
bulence [Q, where the non- universal part of the PDF 
fits with P(Au, r) — r G(^). The same results was 
derived recently with additional assumption by Yakhot 
PH . An interesting point with respect to the possible 
GI breaking mechanisms should be refered. In our cal- 
culation we have shown that C/^ms and the one point 
information has revealed itself because of the matching 
between the inertial range velocity increment PDF and 
the integral scale velocity increment PDF. Because the 
the variance of the velocity increments PDF in the in- 
tegral scales and for the large L is in the order of the 



variance of the one •point PDF we observe such a break- 
down mechanism even with the lack of the forcing con- 
tributions in the Lk m ■ This way of breaking which we 
would like to call as soft GI breaking can be accompa- 
nied with the effects of the forcing where Vrms will be 
entered (JJrms ~ (fc(O)L)^/'^) in the equation of the PDF 
itself. So the second way of the GI breaking which may 
be refered to as a hard GI breaking is linked to the ex- 
plicit dependence of Lkm to fc(0). These parts of the 
Kramers Moyal operator causes some scale dependence 
for the Kramers-Moyal coefficients and the explicit cal- 
culation of the velocity increment PDF in terms of the 
conditional ones is not a trivial task. The intricate part 
of the manipulations is related to the non-commutativity 
of the GI breaking parts of the evolution operator with 
the GI invariant parts so that the scale-ordered expo- 
nential can not be simply operated on the integral scale 
PDF. However it has been shown that these parts would 
cause non-universal behaviours of the amplitudes in the 
velocity increment structure functions. That is for r — > 0, 



5„(r) = A„r^" +3A„_; 



n(n— l)(n — 2) 
' n^B 3 + Cn-3 



-Cn 

(22) 



where Cn — i^n+B) ~ ^ ■ This leads to the non-universality 
of the PDF shapes in the incrtial range jl^. In general 
both of these mechanisms may be relevant for the scaling 
of the inertial range PDF with Au/urms but their clear 
roles would be determined by explicit calculations. 

The A coefficient in eq.(22) is responsible for the scal- 
ing of the structure functions while the B term is an 
infinitesimal coefficient which its value is responsible for 
n independence of the scaling exponents. One should 
note that according to our calculations we have derived 
the following terms just by writing the whole equation 
which is governed over the PDF and then taking all the 
source terms which are proportional to the derivatives of 
the PDF with respect to the angle s. So it is the most 
important result of our calculations which resembles that 
without consulting to the conjectures for introducing the 
scaling terms in the PDF equation all the conjectured 
terms pT[| could be driven just by carefully writing the 
s dependence of PDF and then taking the limit s 1. 
So the A term comes out to be — (rf — 1) = 1 and B 
term approaches to zero as 1 — s^. These are in complete 
coincidence with the proposed values which was derived 
in Yakhot theory for the Burgers turbulence |ll| . Using 
the above equations one can show that the Sn{r) scales 
with r as r^" where C,n = cud — {d — 1) = 1. 

According to Eq.(16) for the universal part in the three 
dimensions we have found — 7/2 while in the nonuni- 
versal part we have got — 3. Comparing the value 
of the density-density exponent ad in universal and non- 
universal regions shows a small deviation in the two parts. 
The dependence of the density density exponent to the 
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specific parts of the velocity increment PDF would be 
revealed if one takes a closer look at the original gener- 
ating function. Considering the eq.(lO), g{r) is referring 
to density-density correlation function conditioned on a 
fixed value of the velocity increment. So it's scaling ex- 
ponent should depend on the behaviour of the velocity 
increment PDF too. 

IV. - DISCUSSION 

The first point which we wish to address is the role 
of viscosity when — s- instead of the inviscid limit. 
It is well known that the viscosity contributions makes 
the generating function equation to become unclosed . 
To our knowledge the only method for dealing with the 
viscosity contribution is the groundbreaking method of 
Polyakov in the one dimensional forced Burgers tur- 
bulence. In this method the overall contribution of the 
viscosity is assumed to generate some other terms con- 
sistent with the symmetries of the original equations and 
their presense just renormalises the terms contributed by 
advective derivative in the equation of generating func- 
tion. Boldyrev generalised the Polyakov's method 
for the Burgers turbulence supplemented with continuity 
equation in one dimension by invoking to the same strat- 
egy from which behaviour of the PDF tails in one dimen- 
sion was given. Keeping the viscous term in the three 
dimensional Burgers equation and with continuity equa- 
tion wc have examined the same closure for the viscous 
contribution. The only relevant term consistent with the 
symmetries developed as a generalisation of the Boldyrev 
work is = a^2 (see for more detail). How- 

ever one can prove that in this framework the strong con- 
straint of homogeniety, i.e. < u >=< u(xi) — u(x2) >= 
will fix the vlaue a = So at least this closure 

has no contribution in PDF equation and the results will 
not be different from the inviscid calculations. 

The second problem which we want to discuss about 
is related to the density probability function in general. 
As far as we know there is one simulation being done by 
Gotoh and Kraichnan jlj] for finding the tail of the den- 
sity PDF in one dimension and for the decaying burgers 
turbulence . For the high density regime they have found 
a power law tail for the density PDF. It is generally ac- 
cepted that in one dimension the mass accumulates in 
the shock regimes and the shock statistics will determine 
the density profile in the stationary regime. However we 
are not aware of such simulations for the forced problem 
even in one dimension. In higher dimensions because 
the nature of the singularities in the velocity flow are 
more involved [3C| ] the simple picture in one dimension 
regarding the mass accumulation in the singularities of 
the velocity profile can not be conducted in a trivial way. 
Apart from the density PDF there is not any simulation 
for investigating any multi-point corralation function of 



the density in higher dimensions. For d = 1 Boldyrev 
JlSt has reported about a simulation on which the expo- 
nent of the two point correlation function of the density 
has predicted to be 2. Any attempt for simulating 
the forced Burgers equation with the continuity equation 
would be valuable for clearing out the outcomes of our 
paper . 

In this paper we have proved the Yakhot conjecture 
about the PDF equation for the multi-dimensional Burg- 
ers turbulence. However he has shown that considering 
B ~ 20 and using the exact result ^3 = 1 as a bound- 
ary condition over final PDF, the incompressible fully 
developed turbulence can be modeled . From his the- 
ory he finds the multiscaling of fully developed turbu- 
lence compatible with the experimental and simulation 
results. In the centrepiece of that work it is emphasised 
that the effect of pressure just renormalises the values of 
B and creates the A term which bahave as source and 
sink. According to our results we have shown that the 
A term could be generated just by considering the angu- 
lar dependences of PDF and even without pressure finite 
scaling exponents are derived and the terms responsible 
for scaling are present in the equation of the stationary 
PDF. It would be illuminating to seek whether, enter- 
ing the effect of pressure just renormalises the values of 
A and B coefficients without adding any new terms in 
the PDF equation. We beleive that the effect of pres- 
sure in the velocity intermittency can be followed if all 
the informations inherent in the correlations of different 
moments of density field are given . Other important 
perspective which we are pursuing is related to the de- 
formation of transvers structure functions and a clear 
picture about the transition between the longitudinal to 
transverse PDF. Consequently this would give a thor- 
ough underestanding about the issue of intermittency in 
the transverse structure functions |^5| . 

Our exact results are also appliciable to the descrip- 
tion of fractal-nature of Interstellar Medium (ISM) if 
we accept that the Burgers equation supplemented with 
continuity eqaution is a good candidate for modelling it 
p2| . In the framework of the Burgers-turbulence the- 
ory of interstellar medium we derive the scaling relation 
M{R) ^ R"^" for the mass on a region of size _R, and the 
value of the dn can be predicted in this framwork which 
is dn = 3/2 . The possible relation of ISM and Burgers 
turbulence will be disscused elsewhere fs^ . 
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